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Abstract 

A theorem on the magnetic energy minimum in a perfect, or ideal, 
conductor is proved. Contrary to conventional wisdom the theorem 
provides a classical explanation of the expulsion of a magnetic field 
from the interior of a conductor that loses its resistivity. It is analo- 
gous to Thomson's theorem which states that static charge distribu- 
tions in conductors are surface charge densities at constant potential 
since these have minimum energy. This theorem is proved here using 
a variational principle. Then an analogous result for the magnetic en- 
ergy of current distributions is proved: magnetic energy is minimized 
when the current distribution is a surface current density with zero 
interior magnetic field. The result agrees with currents in supercon- 
ductors being confined near the surface and indicates that the distinc- 
tion between superconductors and hypothetical perfect conductors is 
artificial. 
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1 Introduction 



In the hterature on superconductivity a distinction between superconduc- 
tors and hypothetical perfect conductors is often made. Even though both 
have zero resistivity, when penetrated by a constant external magnetic field, 
the superconductor expels the field (the experimentally observed Meissner 
effect) but, in the corresponding thought experiment, the perfect conductor 
does not. Here we will prove a theorem on the magnetic energy minimum 
from which one concludes that energy minimization and expulsion of an ex- 
ternal field is the same thing. Therefore, perfect conductors would in fact 
behave like superconductors in their approach to thermodynamic equilib- 
rium. Of course, perfect conductors that are not superconductors exist only 
as thought experiments in textbooks. Nevertheless, textbook thought ex- 
periments should be correct, and unfortunately this is not the case at the 
moment. 

For simplicity we concentrate here on the basic case of type I supercon- 
ductors in a field below the critical field. The outline of this article is as 
follows. We first discuss the historical background that led Meissner and 
subsequent textbook authors astray. We then present our derivation of the 
minimum magnetic energy theorem from a variational principle, after some 
preliminary background on Thomson's theorem and its derivation. Previ- 
ous work on the problem is discussed and then an illuminating example is 
presented in which the energy reduction can be calculated explicitly. Af- 
ter that we briefly discuss the physical mechanism of field expulsion. The 
conclusions are followed by an Appendix which gives further motivation and 
explicit solutions for simple systems. 

2 Meissner effect, diamagnetism, and field ex- 
pulsion 

In the original 1933 work [T] by Meissner and Ochsenfeld0, they stated, with- 
out proof, argument, or reference, that it was understandable that a magnetic 
field would not enter a superconductor, but that the expulsion of a penetrat- 
ing field at the phase transition to superconductivity was unexpected. Since 
then it has been stated innumerable times, in textbooks [21 IH El E] and ped- 

^For a translation of their article into English, see [2]. 
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agogical articles [21 [7], that the expulsion of the magnetic field cannot be 
understood in terms of zero resistivity alone: superconductors are not just 
perfect conductors. The purpose of this article is to show that this is not cor- 
rect so we must first introduce the confused and intricate historical context 
in which the myth arose. 

2.1 The Bohr - van Leeuwen theorem 

In 1911 Niels Bohr [8], and then in 1921, more thoroughly, van Leeuwen [9], 
studied magnetism from the point of view of classical statistical mechanics. 
They reached the conclusion that dia- and paramagnetism could not be ex- 
plained by classical theory. This, so called, Bohr-van Leeuwen theorem is 
often discussed in books on the theory of magnetism, in particular by Van 
Vleck [lOj, but also in more recent texts [11], |12]. The model assumed by 
these authors is a system of classical charged particles, with Coulomb inter- 
actions, in a constant external magnetic field. The key observation is that the 
external magnetic field does not do work on charged particles and therefore 
does not change the energy. Hence no statistical mechanical response. This 
may seem watertight but there are weaknesses. For example, the statistical 
equilibrium distribution is not always determined simply by the energy. If 
there are other conserved phase space quantities, such as angular momentum, 
they may also infiuence the distribution [TH]. Dubrovskii [2] has criticized 
the Bohr-van Leeuwen theorem on these grounds. The main problem though 
is the assumption of a purely external field, and we return to this below. 
Nevertheless, the theorem is still discussed in the literature [151 [13 IHl from 
time to time. 

2.2 The Larmor precession 

There is another, more fundamental, theorem that is in confiict with the 
Bohr-van Leeuwen theorem. Basic Lagrangian dynamics can be used to 
show that a system of particles, with the same charge to mass ratio q/m, 
moving in an axially symmetric external potential, will rotate (precess) with 
the Larmor frequency Q = qB/2'mc, when placed in a weak constant external 
magnetic field parallel to the axis. This is Larmor's theorem [T8l [19] . It is 
easy to show that this precession causes a circulating current which produces 
a magnetic field that screens the external field [2U] . This in accordance with 
Lenz law and is the basis of Langevin's classical derivation of diamagnetism. 
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Feynman et al. [2T] are aware of this discrepancy and indicate that the Bohr- 
van Leeuwen result should not be applied to systems that are free to rotate. 

2.3 The Fock-Darwin Hamiltonian 

The story is, however, even more complicated. A lot of thought has gone into 
understanding the physics behind the Bohr-and Leeuwen theorem. After all, 
the external field makes free charged particles circulate with the cyclotron 
frequency uj = qB/mc, twice the Larmor frequency. Towards the edges of 
the system this should give rise to a net current. It was then realized that 
the forces (or the wall) confining the system has the effect of producing a 
counter current around the edge that cancels this edge effect. The exact 
canceling of the diamagnetic effect by the edge current, however, contradicts 
Larmor's theorem, so the situation is a bit confusing. More detailed studies of 
these effects have been made using the model system of independent charged 
particles confined axially by a parabolic potential in a constant external field 
parallel to the axis. The Hamiltonian for this model is called the Fock-Darwin 
Hamiltonian [221 [23] • In recent years it has been applied to quantum dots. 

One notes that the complete diamagnetism of superconductors is in con- 
flict with the Bohr-van Leeuwen theorem, but not with Larmor's theorem. In 
1938 Welker [23] showed, quantum mechanically, that in a free electron gas, 
confined in a cylinder, the diamagnetic contribution from Larmor rotation 
currents is precisely canceled by the paramagnetic effect on the magnetic 
moments due to orbital angular momenta. This fact led Welker to predict 
an energy gap in superconductors since this would freeze the paramagnetic 
degree of freedom. In spite of all these problems with understanding the 
complete diamagnetism this aspect did not worry Meissner. After all, all 
metals initially exclude an applied magnetic field by appropriate induced 
eddy currents. Due to resistance these subside, but in ideal conductors they 
won't. Only the expulsion of an already present magnetic field at the phase 
transition was therefore considered unexpected by Meissner and Ochsenfeld. 

2.4 The Darwin Hamiltonian 

We note that the findings from the Bohr-van Leeuwen theorem, the Fock- 
Darwin Hamiltonian, and even Larmor's theorem, are not relevant to large 
systems of charged particles. The reason is that considering only an external 
magnetic field is wrong. Unless the charged particles are all at rest they will 
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be sources of magnetic fields as they move. When there are many particles 
this effect, which necessarily is small in few-body systems, easily becomes 
very large ^51 [26], [271 l2H]- For realistic results one must therefore include 
the internal magnetic effects in the relevant dynamics. This was first done 
by Darwin [29]. The Darwin Lagrangian and corresponding Hamiltonian can 
be used to obtain good quantitative understanding of the complete diamag- 
netism, as well as the related effect of rotation called the London moment 

m- 

2.5 Thermodynamics and expulsion 

Already in 1934 Gorter and Casimir ^T] made a thorough analysis of super- 
conductors from the point of view of classical thermodynamics. One result 
of these studies is that the flux expulsion is necessary if the superconducting 
state is to be a thermal equilibrium state. If the final situation depends on 
which order the phase transition and the external magnetic field appeared 
one would have a case of hysteresis, the system would remember its past. 
From this point of view flux expulsion is a natural consequence of approach 
to thermal equilibrium. 

This, of course, does not in itself reveal the actual physical mechanism 
of the expulsion. This is not a subject that has been much discussed in the 
literature with the notable exception of Hirsch. In a number of papers Hirsch 
(see e.g. [321 [331 [311 13S|) speculates on the physical mechanism, and points 
out that the expulsion is not explained by BCS theory. Forrest [2] argues 
that conservation of flux through an ideal current loop, as well as plasma 
physics results on the freezing-in of magnetic field lines, prove that classical 
electromagnetism cannot explain the expulsion. This will be refuted below. 

3 Energy minimum theorems 

Electromagnetic energy can be written in a number of different ways. Here 
we will assume that there are no microscopic dipoles so that distinguishing 
between the D, H and E, B fields is unnecessary. That this is valid when 
treating the Meissner effect in type I superconductors is stressed by Carr 
[36j. Relevant energy expressions are then, 

Ee + E^ = ^f{E' + B')dV (1) 
on J 



5 



l/(,0+l,-.A)dy^ 



(2) 



Iff [ Q{r)Q{r') + ^j(r) ■ j{r') 
2JJ \ \r-r'\ 



) 



(3) 



Here the first form is generally valid while the two following assume quasi 
statics, i.e. essentially negligible radiation. 

According to Thomson's theorem, given a certain number of conductors, 
each one with a given charge, the charges distribute themselves on the con- 
ductor surfaces in order to minimize the electrostatic energy [37J. Even 
though J. J. Thomson did not present a formal mathematical proof for his 
result, the proof of the theorem may be found with great detail in [3S1 [23110] 
and also on its differential form in [H] . A different approach not found in the 
literature, based on the variational principle, is presented here. Thomson's 
result is widely known and has turned to be very useful when applied in sev- 
eral distinct situations. For instance, it was used to determine the induced 
surface density ^2j, and in the tracing and the visualization of curvilinear 
squares field maps [H]. Other applications range from interesting teach- 
ing tools [13] to useful computational methods such as Monte Carlo energy 
minimization [H]. 

For resistive media currents and magnetic fields always dissipate but for 
perfect conductors we will here prove an analogous theorem for the magnetic 
field: The magnetic field energy is minimized by surface current distributions 
such that the magnetic field is zero inside the perfectly conducting bodies, and 
the vector potential is parallel to their surfaces. Since energy conservation is 
assumed, it is clearly only vahd for ideal conductors; otherwise steady current 
exist only if there is an electric field. Previously similar result have appeared 
in the literature [IS], [1^] and we discuss those below. 

3.1 Thomson's theorem 

In equilibrium, the electrostatic energy functional for a system of conductors 
surrounded by vacuum, may be written a^. 



^The ultimate motivation for this specific form of the energy functional, and the corre- 
sponding one in the magnetic case, is that they lead to the most elegant final equations. 
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by combining the electric parts of ([T]) and ([2]) and using E = —Vcj). We 
now split the integration region into the volume of the conductors, V^n, the 
exterior volume, Vont, and the boundary surfaces S, 



dV- f ^(y(PfdV+ f ac^dS, 
JVout on Js 



(5) 



where a is the surface charge distribution. 

We now use, (V0)^ = V ■ (0V0) — 4>V^4>, and rewrite the divergencies 
using Gauss theorem. The energy functional then becomes. 



+ 


1 




dV 




(f)dV 










JVout Sit 






(j(f) 


1 


pn ■ 


(v+0 - V-(f)) 


dS, 













(6) 



where V"*" and V~ are the gradient operators at the surface in the outer and 
inner limits, respectively. Since the total charge in the conductor is limited 
and constrained to the conductor volume and surface, one must include this 
restriction by introducing a Lagrange multiplier A. Therefore the infinitesi- 
mal variation of the energy is. 



SE, 



+ 



6(j)[g + —V'(j)]+6g{(j)-X) 
a + ^n- (V" 



dV + 



An 



-V> dV 



A)| d^. 



From this energy minimization, the Euler-Lagrange equations become: 



= -Ang 
(h = \ 



(7) 



(8) 



S 



= A 



out 







(9) 
(10) 



According to these equations the potential is constant inside the conductor 
in the minimum energy state. This means the electric charge distribution 
is on the surface and in an equipotential configuration. This concludes the 
proof of Thomson's theorem. 
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3.2 Minimum magnetic energy theorem 



A similar procedure will now be applied to the magnetic field. We write the 
magnetic energy functional for a time independent magnetic field is written 

as, 

i?^ = / [- j • A - ^ (V X Af] dV, (11) 

i.e. as two times the form ([2]) of the magnetic energy minus the form ([T]), 
using S = V X A. As before we split the volume into the volume interior 
to conductors, the exterior vacuum, and the surface at the interfaces, and 
write. 



1 



1 



= I -J . A - (V X Af dV 

JVin LC OTT 

-/ —{VxAYdV+-[k-AdS, 

Jvout 8n c Js 

where k is the surface current density. We now use the identity, 
(V X A)^ = V ■ [A X (V X A)] + A ■ [V X (V X A)] , 



(12) 



(13) 



and then use Gauss theorem to rewrite the divergence terms. The energy 
functional then becomes: 



f |lj.A-^A.[Vx(VxA)]jd\/ 

- / —A - [V X (V X A)]dV 

JVout Stt 



(14) 



+ 



-k-A-—A 

sic Stt 



nxfV+xA-V xA 



dS. 



As in the electric field case, some constraints should be imposed. In particu- 
lar, due to charge conservation, the electric current must obey the continuity 
equation both in volume and at the surface |47j : 



Vs-k 








(15) 
(16) 



where V5 is the surface gradient operator. Notice that the constraints are 
local, not global. In other words, the Lagrange multiplier is not constant 
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but a space function A(r). Therefore, infinitesimal variation of the magnetic 
energy gives, 

c5E^ = j^^ |(5A . (^j - ^ [V X (V X A)]) + 5j ■ (A - VA)| 

- / 5A ■ ^ [V X (V X A)] dV (17) 

JVout 47r 

+ ^ I^A ■ A; + X (V+ X A - V" X a) + 5fc • (A - V5A)| d5. 

Equating this to zero we find that, 

r Vx(Vx A) = VxS = ^j, 
A = VA , 



and, 



and. 



/ fc = t^x(V+xA-V-x A) , 

1 A = V5A, ^'^^ 



Kut V X (V X A) = V X S = 0, (20) 



are the Euler-Lagrange equations for this energy functionaL 

According to eq. (fTSjl B = V x VA = 0, so the magnetic field must be 
zero in V^^. Consequently also the volume current density is zero, j = 0, 
inside the conductor, in the minimum energy state. 

3.2.1 Surface currents 

Now consider the results for the surface, eq. (fT9l) . Our results from Vin show 
that X A = 0, so the equation reads, 

k = —h X (V+ X V5A) . (21) 

Let us introduce a local Cartesian coordinate system with origin on the sur- 
face, such that the surface is spanned by x, y with unit normal fi = z = x xy. 
Assuming that the surface is approximately flat we then have that, 

^s = x— + y—, and, = x— + y— + n—- = V s + n—-. (22) 
ox oy ox oy oz^ oz^ 
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Since, A = V^A, the vector potential is tangent to the conducting surface 
and we get, 



47r 

— k = n X 

c 



Vs + n- 



d 



n X 



n 



d 
9i+ 



X X 



dz+ 

dX ^ ^dX\ 
dx ^ dy J 



X V5A(x,y,2;; 
. ( . d 



n X \ n 



dz+ 



X A 



(23) 
(24) 



for the surface current density. Rewriting the triple vector product we find, 



d 



c dz 



[(n ■ A)n — {n ■ n)A\ 



dA 

"dz+' 



(25) 



so the surface current density is parallel to the outside normal derivative of 
the vector potential. We note that this agree with the well known result ^0] , 



k = h X (B, 



B-) 



(26) 



for the case of zero interior field {B_ = 0). 



3.2.2 External fields 

Even though the proof of the theorem does not include the presence of an 
external magnetic field, the result can be generalized. Assume that the mag- 
netic field is produced by currents at infinity. These currents demand a 
non-zero surface integral at infinity in eq. f|T^ which does not affect the fi- 
nal result, however. Alternatively one can include large perfectly conducting 
Helmholtz coils (tori) in the system and let the rest of the system be small 
and situated in the constant field region of the coils. In section H] we instead 
treat this case by means of an illuminating example. 



3.3 Previous work 

The fact that there is an energy minimum theorem for the magnetic energy of 
ideal, or perfect, conductors, analogous to Thomson's theorem, is not entirely 
new. In an interesting, but difficult and ignored, article by Karlsson [35] such 
a theorem is stated. Karlsson, however, restricts his theorem to conductors 
with holes in them. In the electrostatic case charge conservation prevents the 
energy minimum from being the trivial zero field solution. In the magnetic 
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ideal conductor corresponding conservation law is the conservation 

of magnetic flux through a hole |1H]. But, as long as one conductor of the 
system has a hole with conserved flux there will be a non-trivial magnetic 
field. To require that all conductors of the system have holes, as Karlsson 
does, seems to us unnecessarily restrictive. One of his main results is that a 
the current distribution on a superconducting torus minimizes the magnetic 
energy. 

A result by Badfa-Majos [46j comes even closer to our own and we outline 
it here. The current density is assumed to be of the form, 

j = qnv, (27) 

where q is the charge of the charge carriers and n is their number density. 
The time derivative is then given by, 

dj an dv an f ^ a \ q'^n „ o . „ 

-± = :^m— = — [qE + -vxB]=—E + —jxB, (28) 

dt m dt m \ c J m mc 

assuming that only the Lorentz force acts (ideal conductor). We now recall 
Poynting's theorem |l9j for the time derivative of the field energy density of 
a system of charged particles, 

d fE^ + B^\ . _ „ / c 



j -E-V ■ [ —E X B). (29) 



dt \ 8iT J KAtc 

The first term on the right hand side normally represents resistive energy 
loss. Here we use the result for E from eq. (l28l) . 



and get. 



B, (30) 

g^ndt qnc 



. ^ m dj . d / m , , 

j.E=—/-j = — [-^j']. (31) 
g^n dt dt yZq'^n I 



This is thus the natural form for this term for perfect conductors. We insert 
it into (129|) . neglect radiation, and assume that E"^ -C B^. This gives us. 



1 „o m 




^ ikn ^ °- '''' 
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Finally inserting, j = (c/47r)V x B, here, gives. 



E 



B 



1 

Stt 



V JVin 47rg^n 



(V X B)'^dV 



(33) 



for the conserved energy, after integration over space and time. 

Badfa-Majos [46] then notes that this energy functional implies flux ex- 
pulsion from superconductors. Variation of the functional gives the London 
equation [19] . 

1 TflC^ 

B + X (V X = 0. (34) 

Badfa-Majos, however, does not point out that this classical derivation of flux 
expulsion is in conflict with text book statements to the effect that no such 
classical result exists. Further work by Badfa-Majos et al. [50] on variational 
principles for electromagnetism in conducting materials should be noted. 



4 Ideally conducting sphere in external field 

We now know that magnetic energy minimum occurs when current flows only 
on the surface and the magnetic field is zero inside. Let us consider a perfectly 
conducting sphere in a fixed constant external magnetic field. Here we will 
calculate the surface current needed to exclude the magnetic field from the 
interior and how much the energy is then reduced. In order to exclude a 
constant external field B from its interior the currents on the sphere must 
obviously produce an interior magnetic field —B, thereby making the total 
field zero in the interior. It is well known that a constant field is produced 
inside a sphere by a current distribution represented by the rigid rotation of 
a constant surface charge density. 

4.1 Energy of the external field 

For the total magnetic field to have a finite energy we can not assume that the 
constant external field extends to infinity. Instead of using Helmholtz coils 
we simplify the mathematics and produce our external field by a spherical 
shell of current that is equivalent to a rigidly rotating current distribution on 
the surface. This can be done in practice by having as set of rings represent- 
ing closely spaced longitudes on a globe with the right amount of currents 
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maintained in each of them. Such a spherical shell of rigidly rotating charge 
produces a magnetic field that is constant inside the sphere and a pure dipole 
field outside the sphere (see Fig. [1]): 




B{r) 



2m 

3(m ■ f)f — m 



for r < R 



for r > R 



(35) 



Here r = |r| and the center of the sphere is at the origin. If Q is the total 
rotating surface charge and w its angular velocity, 

QR' 



m 



3c 



(36) 



see eq. below. It is now easy to calculate the magnetic energy of this 
field. One findfl 



En = — 



B'dV 



r<R 



r>R 



1\ m 



m 



Inside this sphere, which is assumed to maintain a constant magnetic field 
in its interior, we now place a smaller perfectly conducting sphere. 



■^If (|36|) is inserted for m here we get Eq = [uj / c)'^ RQ'^ / 9 which is equal to of Eq. 
([56|) below, for ^ = 1, corresponding to surface current only, as it should. 



13 



4.2 Magnetic energy of the two sphere system 

We assume that the small sphere in the middle of the big one has radius a < R 
and that it also produces a magnetic field by a rigidly rotating charged shell 
on its surface. We denote its dipole moment by so that its total energy 
would be, 

far from all other fields, according to our previous result fl57|) . We now place 
the small sphere inside the large one and assume that makes an angle a 
with m, 

m ■ rUa = mniaCosa. (39) 
The total energy of the system is now, 

E = ^ f{B + B^fdV = Eo + E, + E„ (40) 

where the interaction energy is, 

E, = ^Jb- BadV (41) 

The integral here must be split into the three radial regions: < r < a, 
a < r < R, and R < r. The calculations are elementary using spherical 
coordinates. The contribution from the inner region is 

En = ^f B.BAV = ^'^cosa (42) 

47r Jr<a 6 R^ 

The middle region where there is a dipole field from the small sphere in the 
constant field from the big one contributes zero: Ei2 = 0. The outer region 
gives Ei3 = {2 / ?))mm a CO?, a / R? . Summing up one finds, 

Ei = 2—^ cos a, 43 
R-^ 

for the magnetic interaction energy of the two spheres. 

4.3 Minimizing the total magnetic energy 

The total magnetic energy of the system discussed above is thus, 

2 2 

Eo + Ea + E, = E{ma, a) = — + ^ + 2-^ cosa. (44) 
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We assume that m and are positive quantities. This means that as a 
function of a this quantity is guarantied to have its minimum when cos a = 
— 1, i.e. for a = 7T. Thus, at minimum, the dipole of the inner sphere has the 
opposite direction to that of the constant external field. 

Now assuming a = vr we can look for the minimum as a function of rria- 
Elementary algebra shows that this minimum is attained for. 



ma = y—j m = mmin- (45) 

The magnetic field in the interior of the inner sphere is then, 

_ _ 2m 2ma /2m 2mmin\ m , , 

B + S, = — + — p _ = 0, 46 



so it has been expelled. The minimized energy (jH]) of the system is found to 
be, 

2 

771 

-E'min = E{r7l^i^, 7l) = — 



a \3 

r) 



(47) 



The relative energy reduction is thus given by the volume ratio of the two 
spheres. Using B = 2m/ R^ from (135|1 the energy lowering can also be 
expressed as follows: 

^^-^--R.R.---'-\^)\-Y-)- 

This result is independent of the radius of the big sphere introduced to pro- 
duce the constant external field. It shows that the energy lowering corre- 
sponds to three times the external magnetic energy in the volume 47ra^/3 of 
the perfectly conducting sphere. 



5 The mechanism of flux expulsion 

Assume that a resistive metal sphere is penetrated by a constant magnetic 
field. Lower the temperature until the resistance vanishes. How does the 
metal sphere expel the magnetic field, or equivalently, how does it produce 
surface currents that screen the external field? According to Forrest [2] this 
can not be understood from the point of view of classical electrodynamics 
since in a perfectly conducting medium the field lines must be frozen-in. This 
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claim is motivated thus: When the resistivity is zero there can be no electric 
field according to Ohm's law, since this law then predicts infinite current. 
But if the electric field is zero the Maxwell equation, V x — dB/dt = 0, 
requires that the time derivative of the magnetic field is zero. Hence it must 
be constant. 

Ohm's law is, however, hardly applicable in this situation. The system of 
charged particles undergoes thermal fluctuations and these produce electric 
and magnetic fields. The electric fields accelerate charges according to the 
Lorentz force law. In the normal situation these currents stay microscopic. 
When there is an external magnetic field present the overall energy is lowered 
if these microscopic currents correlate and grow to exclude the external field. 
According to standard statistical mechanics the system will eventually relax 
to the energy minimum state consistent with constraints. We note that 
Alfven and Falthammar |51] state that "in low density plasmas the concept 
of frozen-in lines of force is questionable" . 

Another argument by Forrest [2] is that the magnetic fiux through a per- 
fectly conducting current loop is conserved. Since one can imagine arbitrary 
current loops in the metal that just lost its resistivity with a magnetic field 
inside, the field must remain fixed, it seems. Is is indeed correct that the 
fiux through ideal current loops is conserved, but the actual physical current 
loop will not remain intact unless constrained by non-electromagnetic forces. 
On a loop of current that encloses a magnetic fiux there will be forces that 
expand it [ID]. This is the well know mechanism behind the rail gun, see e.g. 
Essen |52]. All the little current loops in the metal will thus expand until 
they come to the surface where the expansion stops. In this way the interior 
field is thinned out and current concentrates near the surface. So, the fiux is 
conserved through the loops, but the loops expand. 

6 Conclusions 

We have now reached the conclusion that the Meissner effect, both the perfect 
diamagnetism and the flux expulsion, are consequences of classical electrody- 
namics, together with reasonable assumptions about the system of charged 
particles. The reason that this has not been recognized before is a widespread 
ignorance of the properties and effects of the magnetic interaction energy 
between moving charged particles. The majority of the physics community 
tends to think of magnetism as either due to microscopic dipoles in a medium 
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or as due to an imposed external magnetic field. That the translational mo- 
tion of the charged particles of the system also produces magnetic fields, 
which in turn change their motion, is normally neglected. When this can not 
be neglected one must find a self consistent simultaneous solution for cur- 
rents and magnetic fields [53] . More awareness of these facts will hopefully 
clear up a lot of confusion in the future. 



A Appendix 

A.l Magnetic energy minimization in simple one de- 
gree of freedom model systems 

We now define simple one degree of freedom model systems and use them 
to illustrate how the minimum magnetic energy theorem works. We take 
systems in which all energies can be calculated exactly so that energy mini- 
mization amounts to minimizing a function of a single variable. The systems 
are both related to a system used by Brito and Fiolhais [13] to study electric 
energy. 



A. 1.1 Magnetic energy of coaxial cable 

The minimum magnetic energy theorem can be illustrated in such a simple 
system as a coaxial cable. The cable can be modelled by an outer cylindrical 
conducting shell with radius b, carrying an electric current J, and a concentric 
solid cylindrical conductor with radius a < b, carrying the same electric 
current in the opposite direction. We now divide the inner current into a 
surface current Is = {1 — a.t p = a and a bulk current 1^ = 1]I in the 
inner volume, < p < a, so that I = 1^ + ly. We use cylindrical coordinates, 
p, (p, z, and obtain the magnetic field. 



B 
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if a < p < b 



(49) 



b < p 
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using Ampere's law. Thus, the magnetic field energy for a length L of the 
cable is given by: 



Svr Jv 



2IV 
vr V c / 



b 



LP 



4 \a 



(50) 



This magnetic energy reaches its minimum for zero bulk current, 77 = 0, 
corresponding to surface current only and zero field for < p < a. 



A. 1.2 Current in sphere due to rigidly rotating charge 

Consider an ideally conducting sphere of radius R. Assume that there is a 
circulating current in the sphere which can be seen as the rigid rotation of a 
charge Q evenly distributed in the thick spherical shell between r = a < R 
and r = R. The charge density. 



r 



SQ 



for < r < a 



for a < r < R 



for i? < r 



(51) 



is assumed to rotate with angular velocity uj = uo z relative to a an identical 
charge density of opposite sign at rest. The current density is then. 



Q{r) (JL> X r, 



(52) 



and the current, / = ^t^, passes through a half plane with the z- axis as 
edge. 

The vector potential produced by this current density can be found using 
the methods of Essen ^20j, see also SO]. If we introduce ^ = a/R we 
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find, 



A{r) = —{u}xr)-i 



for < r < a 



/ 5 

.^/. ..V 7^-37^-2^ for a<r<i? (53) 
10(1 - ^3) I i?3 R5 J - - ^ ^ 



(1 - 

[ 5(1 - r3 



for < r 



The parameter ,^ = a/R is zero, ^ = 0, for a homogeneous ball of rotating 
charge, while ^ = 1 corresponds to a rotating shell of surface charge. Com- 
paring with the vector potential for a constant field, A — ^Bq x r we see 
that, 

gc^ (i-a _Q^ (1 + 

° ci?(i-a ci?(i+e+a ^ ^ 

is the field in the central current free region <r < a. 



A. 1.3 Magnetic energy of rotating spherical shell current 

We now calculate the magnetic energy of this system using the formula. 



Performing the integration using spherical coordinates gives. 



where. 



2 + 4^ + 6^2 ^ + 10^^ + 5^^ 



(55) 



(56) 



(57) 



35(1 

is a function of the dimensionless parameter ^. 

This expression for the energy is the Lagrangian form of a kinetic energy 
which depends on the generalized velocity cu — (p, 



(58) 



19 



Since the generalized coordinate if does not appear in the Lagrangian 
the corresponding generahzed momentum (the angular momentum), 



(59) 



is a conserved quantity. The corresponding Hamiltonian form for the energy 
is then Hm = Pif<p — L^, expressed in terms of p<^, 



(60) 



The function l/f{C,) is plotted in Fig. [2l We now consider the two energy 




0,0 0,2 0,4 0,6 0,8 1,0 

Figure 2: A graph of the function l//(^) which is proportional to the Hamiltonian 
form of the magnetic energy (j60p of our model system. Note that ^ = corresponds 
to volume (bulk) current and ^ = 1 to pure surface current. 

expressions ( I58l) and (1601) separately. 

Case of constant current: We first consider the case that the current, 
I = Quj/2tt, is constant. Changing ^ then means changing the conductor 
geometry while keeping a constant total current, or, equivalently, angular 



20 



velocity uo = fp. One might regard the total current as flowing in a con- 
tinuum of circular wires. Changing ^ from zero to one means changing the 
distribution of these circular wires from a bulk distribution in the sphere to 
a pure surface distribution, while maintaining constant current. According 
to a result by Greiner [51] a system will tend to maximize its magnetic en- 
ergy when the conductor geometry changes while currents are kept constant. 
This has also been discussed in Essen ■ In conclusion, if currents are kept 
constant the magnetic energy (158!) will tend (thermodynamically) to a stable 
equilibrium with at a maximum value and we note that this corresponds to 
a pure surface current ^ = 1. 

Case of constant angular momentum: Assume now that we pass 
to the Hamiltonian (canonical) formalism. Thermodynamically this type of 
system should tend to minimize its phase space energy fl60|) in accordance 
with ordinary Maxwell-Boltzmann statistical mechanics. As a function of ^ 
this Hamiltonian form of the energy clearly has a minimum at ^ = 1, see 
Fig. [2], corresponding to pure surface current. In this case therefore there 
will be current density only on the surface in the energy minimizing state. 
This is in accordance with the our minimum magnetic energy theorem. It 
is notable that both the assumption of constant current and the assumption 
of constant angular momentum lead to a pure surface current density as the 
stable equilibrium. 

A. 2 Explicit solutions with minimum magnetic energy 

Here we present two explicit solutions for the current distributions and mag- 
netic fields that minimize the magnetic energy. The solution for a torus has 
been found several times, probably first by Fock [53], but, independently, 
several times since then, see e.g. [SH E?! EHl EH EDI EI], so we do not re- 
peat it here. Karlsson ^5], however, was probably the first to notice that 
the solution minimizes magnetic energy for constant flux. Dolecek and de 
Launay [IB] verified experimentally that a type I superconducting torus be- 
haves exactly as the corresponding classical perfectly diamagnetic system for 
field strength below the critical field. Here we treat two simpler cases involv- 
ing constant external fields: a cylinder in a transverse magnetic field, and a 
sphere. 
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A. 2.1 Cylinder in external perpendicular magnetic field 

Consider an infinite cylindrical ideal conductor with radius R in a external 
constant perpendicular magnetic field. To get the vector potential one must 
solve the following differential equation, 

V X S = V X (V X A) = 0. (61) 

To solve this one should look for the symmetries of the system. We assume 
that the external constant magnetic field points in the ^/-direction and that 
the cylinder axis coincides with the z-axis. There will then be no dependence 
on the ^-coordinate so the magnetic field is, 

Moreover, due to the symmetry of the system, the z-component of the mag- 
netic field must be zero, 

|(p^.)-|f-0. (63) 



These assumptions and constraints transform eq. (161]) into, 

V X S - - {—^^ ^ - f64l 

\p2 df^"^ dp^ p dp J 

which simply is Laplace equation in cylindrical coordinates {p,(f,z). Before 
writing down the general solution, let's analyze the boundary conditions. As 
p — > oo, the magnetic field must approach the external one: = Bq]) = 
i?o(<^cos(/9 + p sin if). Furthermore, since the magnetic field is zero inside 
the perfect conductor, one concludes from eq. fl62|) that the vector potential 
vector must be constant inside the cylinder. Therefore, the solution is, 

Az = const. -|- Bq ^ pj cos if, (65) 

for p > R. The magnetic field outside the cylinder becomes, 

B = p Bq (l - ^\ sin(p + (p Bq (l + cos (f, (66) 
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which imphes that, 

B{p = R) = 2Bo cos cp ip (67) 

The magnetic field on the cyhnder's surface determines the surface current 
according to eq. (1261) . so we get, 

c 

k = — i?ocosv9 z. (68) 

The total current obtained through integration of the surface current is zero 
as expected, otherwise the energy would diverge. A more detailed analysis 
on this problem can be found in [22] • 



A. 2. 2 Superconducting sphere in constant magnetic field 

Similar calculations can be performed for a superconducting sphere with 
radius i? in a constant external magnetic field pointing in the direction of 
the 2;-axis. 

Similarly to the cylinder case, eq. ( l6Ti) is simplified a lot due to the sym- 
metries of the system. Since the external constant magnetic field points along 
the z-axis, there won't be any dependencies on the (p coordinate and the mag- 
netic field along this coordinate must be zero. Therefore, the magnetic field 
simplifies to, 

B = {A^ sin e)f---^ (V) 9, (69) 

r sm u GO r or 

where we use spherical coordinates r, 6, ip. Again, using these assumptions 
and constraints eq. fl6Tl) becomes. 



V X B = -- 

r 



d ( d \ 1 d ( 1 d . 
Q-r Yr ^"^^^ j + ; M 9^ 



^ = (70) 



Since the magnetic field is zero inside the sphere, eq. (l69i) implies that the 
vector potential has the form, 

A^{r <R) = (71) 
r sm fc' 

where C is a constant. To prevent the vector potential from diverging at 
r = and 6* = 0, the constant C must be zero. Furthermore, as r — oo, the 
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magnetic field must go to the external field, Bq = Bqz = Boi^f cos 6* — ^ sin ^) . 
Therefore, the solution of eq. (17U]) for this case is, 

A^{r>R) = ^(^r-^ym9, (72) 

which leads to the following magnetic field outside the sphere, 

^ cos^-^^fio 1 + ^- 



S = f5o ( 1 - — cos^ - ^fio 1 + t:— sin^. (73) 



The magnetic field at the sphere surface is thus. 



B = -^BQ^ind d. (74) 



One notes that this is the same field as that of section 14.21 at the surface of 
the inner sphere when energy is minimized. 

Using eq. fl26l) . the surface current density becomes, 

fc = — ^5osin0(^ (75) 

Unlike the infinite cylinder in a perpendicular external field, the sphere must 
have a total non-zero electric current, / = ^RBq, to keep the magnetic field 
from entering. A similar approach to this problem can be found in [63] . 
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